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Abstract 

We study an intrinsic curvature model defined on fixed-connectivity triangulated lattices 
enclosing a spherical core by using the canonical Monte Carlo simulation technique. We 
find that the model undergoes a discontinuous transition of shape transformation between 
the smooth state and a collapsed state even when the core radius R is sufficiently large; 
the transition depends on R. The origin of the multitude of transitions is considered to be 
a degeneracy of the collapsed states. We also find that the Gaussian bond potential S \ IN, 
which is the sum of bond length squares, discontinuously changes at the transition. The 
discontinuity in S \IN implies a possibility of large fluctuations of the distance between 
lipids, or the density of lipids, in biological membranes such as giant vesicles or liposomes 
enclosing some materials. 
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1 Introduction 



A lot of attention has been paid to the phase structure of surface models for strings 
and membranes nil2l3l4l5| . The phase transitions of the models are closely con- 
nected to a variety of morphology in membranes [7 8.9I10I11I12L The Hamilto- 
nian is defined by the squared mean curvature or the extrinsic curvature on two- 
dimensional smooth surfaces embedded in [I13I14II . and the statistical mechan- 
ical models are obtained by discretizing the Hamiltonian on the triangulated sur- 
faces II15I16I17I18I19I . In numerical studies, phantom surfaces are always assumed 
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as membranes for computational efficiency; a surface which is allowed to self- 
intersect is called the phantom surface. The phase structure of the phantom surfaces 
has been extensively studied and is partly clarified such that the spherical surface 
models undergo a first-order transition of shape fluctuations [12012 1 122I23I24I25 (26ll . 

Membranes are always understood as boundary surfaces separating interior materi- 
als from exterior environment. Membranes can shrink/swell if the interior material 
is liquid. Liquid can leak out/in through pores or by some other reasons, and con- 
sequently the membrane fluctuates as a two-dimensional surface. Thus, enclosed 
materials, such as liquid, and their property are considered to be included in the 
surface model as the bending rigidity. To the contrary, hard materials can not sim- 
ply be included in the surface model. Surface fluctuations of membranes are con- 
siderably suppressed by such hard materials. Membranes which enclose the hard 
materials never shrink. 

Thus, the phase structure of membranes enclosing hard materials is expected to be 
different from the phase structure of membranes without the hard materials. It is 
possible that the collapse phenomenon does not occur on the surfaces enclosing a 
spherical core. The surface is prohibited from changing the orientation upside down 
owing to the core inside the surface even when the surface is phantom. 

Therefore, it is interesting to study the phase structure of surface models with a 
spherical core and to see the dependence of the phase structure on the core radius 
R. In this paper, we study a fixed-connectivity spherical surface model in [|27ll with a 
spherical core by using the canonical Monte Carlo (MC) simulation technique. The 
model is an intrinsic curvature model [|28|29I30II and is different from the extrinsic 
curvature ones in [2 0.21.22i23l24l25l26L The intrinsic curvature is independent 
of the surface shape, and hence its value is expected to be almost independent of 
whether the surface encloses the spherical core or not. 



2 Model 

2. 1 Fixed-connectivity phantom surface models 

In this subsection, we make a brief survey of the triangulated surface models and 
their phase structures. 

The Hamiltonian S is given by a linear combination of the bond potential S i and 
the curvature energy S 2 such that S = S i +bS 2, where b[kT] is the bending rigidity. 
The continuous Hamiltonian is discretized on triangulated surfaces of sphere topol- 
ogy in for numerical studies. Muenkel and Heermann reported in Ref. [20J that 
the model undergoes a phase transition between the smooth phase and the crumpled 
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phase on both of the sphere and a torus in R^, where the Hamihonian S includes an 
intrinsic curvature term, and moreover that a self-avoiding surface model undergoes 
a discontinuous transition [21 J. The role of 5 1 is to make the mean bond length con- 
stant. Therefore, a Lennard- Jones type potential can be assumed as S i [21 22J and, 
a hard wall potential is also possible [|26l| . In Ref. lfTSl of Kantor and Nelson, a hard 
wall and hard core potential is assumed as 5 1; the model is well known as the ball- 
spring model. Gompper and KroU extensively studied the surface model including 
the ball-spring model on both phantom and self-avoiding surfaces [ 161171181 . On 
the other hand, the role of 5 2 is to make the surface smooth at sufficiently large 
finite b. 

In the model of Helfrich and Polyakov, S 1 is given by the Gaussian bond potential; 
•^1 = Tjij(Xi-Xj)^, where X,(e R^) is the vertex position, and ^2 = Zz/l^w, • 
n^), where n, is a unit normal vector of the triangle i. We call this model as the 
canonical surface model. The canonical model, as mentioned in the Introduction, 
was reported to have a first-order phase transition separating the collapsed phase 
from the smooth phase [|25l . Note also that the canonical surface model is expected 
to have a continuous transition from theoretical calculations on the basis of the 
renormalization group [3.4.5.6,1 . 

The intrinsic curvature model is obtained by replacing S 2 with the intrinsic curva- 
ture energy 52"' = - 2, log(5,/2;r), where is the sum over angles of triangles 
meeting at the vertex /. This term corresponds to a weight 11,^°" = 

exp(cr2, log^,) for the integration of the partition function, where ^, is the co- 
ordination number of the vertex / and cr = d/2 = 3/2 [[351 . As mentioned above, 
surface models defined by a Hamiltonian which includes the term -(3/2) 2; log 
are extensively studied H20I28I29I30II . while the constant parameter cr is extended 
to a variable one A, and a model with such an intrinsic curvature term A 2; 
was also extensively studied [28.1 . 5'2"' = Xi l^i-^M or S2^ = Yji(Si-'2-n)^ can also be 
assumed as the intrinsic curvature energy, where the deficit angle 6i-2n is exactly 
zero on the planar lattices. It is apparent that 5, = In on the cylindrical lattices if the 
triangulation is sufficiently fine. Thus, 5 2"' does not always depend on the surface 
shape in contrast to the extrinsic curvature energy. However, numerical studies in- 
dicate that the intrinsic curvature model undergoes a first-order transition of shape 
transformation; not only the collapsed phase but also the smooth phase is stable 
[|77]| . This is a non-trivial result in the sense that S^2^ uniquely determines the sur- 
face shape even though 5, is not always dependent on the surface shape. 

The phase transition of the intrinsic curvature model can be seen on spheres, tori, 
and disks Il27l31|32i Hence the transition is independent of these three types of 
topology. 

As mentioned above, 5 2"' is zero not only on the planar surface but also on the tubu- 
lar surface. Thus, the smooth phase is considered to be degenerate in the model on 
the disk surface, because the planar disk can turn a tubular surface without a cost 
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of 5* 2"'. However, the smooth surface is almost planar on the disk surface ||32]| . This 
implies that the smooth phase is "almost degenerate". This allows us to consider 
that the collapsed phase is also degenerate or almost degenerate. Many of the col- 
lapsed surfaces are expected to belong to a single value of 5 or to a narrow and 
continuous range of 5 2"'. We numerically check this expectation in this paper and 
show that it is correct. 

2.2 Model with a hard-core potential 




(a) (6) 

Fig. 1. (a) A schematic view of a surface enclosing a spherical core of radius R, (b) 
the smooth surface configuration of triangulated sphere of radius RoiN), which satisfies 
S i/N = 3/2. The surface in (b) is expected to appear in the limit of a— > 00. 

In this subsection, we define an intrinsic curvature model with a hard core potential. 
Figure [TJa) shows a spherical surface enclosing a spherical core of radius R. The 
surface size is always greater than the core size even in the collapsed phase. 

The triangulated spherical lattice in this paper is identical to that in [|25l and is 
constructed from the icosahedron by splitting the edges (faces) into £ {f-) pieces. 
The total number of vertices is thus given by = l0€^+2, the total number of bonds 
Nb and the total number of triangles Nj are given by Nb = 30{^ and Nj = 20i^. The 
coordination number ^ is ^ = 6 almost everywhere on the lattice excluding the 12 
vertices, which are the vertices of the icosahedron and of ^ = 5. 

The partition function Z is defined by 

Z= ]~[jX,exp[-5(X)], (1) 

^ i=i 

where J Y\f=i dXi is the BA'^-dimensional integrations, and X, R^^ is the vertex 
position /. The Hamiltonian S{X) denotes that S depends on the surface shape 
X{= {Xi,X2, ■ ■ ■ ,Xn]). The Hamiltonian is defined by a linear combination of the 
Gaussian bond potential 5 1, the deficit angle energy 5 2"', and the hard-core poten- 
tial U{R) such that 



4 



S(X) = S I + aSf + U(R), 

S:=J](Xi- Xjf , S'f = -J] log {6i/2n) , (2) 

('7) i 

U(R) = Ui{R). 

i 

Yj(ij) in S 1 denotes the sum over bonds (ij). The symbol a[kT] is the curvature 
coefficient, and 5, in 5 2"' is the sum over angles of triangles meeting at the vertex 
/; 6i was introduced in the previous subsection l2.1[ 2, in S^2^ denotes the sum over 
vertices i, and Ui(R) in U (R) is given by 

Ui(R) = ' (3) 
I 00 (otherwise) . 

We should note that the center of the core is fixed at the origin of R^, while the 
center of mass of the surface is not fixed in the integrations of Z in contrast to the 
integrations of Z of the models without the core. The translational zero mode of the 
surface can influence the results if the surface configurations are partly dominated 
by such configurations that are identical to each other up to translation. However, 
such configurations are not expected, moreover, the core protects the surface from 
the translation. The collapsed surface is at least prohibited from the translation, 
because the size of collapsed surface is not so larger than the core size. As for the 
smooth surface, the surface translation is also expected to be almost suppressed. 

We should note that the scale invariance of the partition function is broken by the 
term U(R) [12J. As a consequence, the relation (S JN) = 3/2, which holds in the 
model without the cores, is expected to be violated because of the term U(R) and 
the scale invariance of the partition function [.12.1 . 

Table 1 

The assumed values of N, Rq, and R. The core size R ranges from R(50%) to R{20%). 



N 


^0 


R{50%) 


7?(30%) 


R{20%) 


1442 


7.046 


3.523 


2.1138 


1.4092 


2562 


9.39 


4.695 


2.817 


1.878 


4842 


12.9 


6.45 


3.87 


2.58 


6252 


14.66 


7.33 






8412 


17.01 


8.505 






10242 


18.76 


9.38 






12252 


20.52 


10.26 







We assume several diff"erent values of R such that R/Rq = 50%, R/Rq = 30%, R/Rq = 
20%, and R/Rq = 0%. The symbol Rq denotes the radius of smooth sphere, where 
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the bond length squares {Xi-Xjf satisfies Y.(ij){Xi-Xjf IN(= S ^N) = 1.5. Rq is 
identical with the radius of the equilibrium surface in the limit of a ^ oo and is 
determined only by A^; Rq is a parameter that depends only on A'^. Figure[2b) shows 
a sphere section of radius Rq. The core radius R is also determined by both and 
the ratio RIRq. The assumed values of Rq and R are shown in Table [TJ Note also 
that Rl/N, and hence R^/N also, becomes a constant independent of N including 
the limit of £». In fact, the area 4-7tRI of the surface of radius Ro is proportional 
to A^, because the mean bond length squares is a constant, which is independent of 
A^ due to the fact that S i /A^(= 1.5). 



3 Monte Carlo technique 

The vertex position X on the triangulated surface is sequentially updated by using 
the canonical Metropolis MC technique. A new position X' = X + 6X is accepted 
with the probability Min[l, exp(-(55)], where 6S = 5 (new) - 5 (old). The small 
variation 6X is randomly generated in a sphere of a given radius, which is fixed to 
maintain approximately 50% acceptance rate of the Metropolis accept/reject. The 
Metropolis algorithm is applied only to vertices, of which X' is located outside the 
core. New positions X' inside the core are rejected. This acceptance rate of X' is 
almost independent of R and, almost all X' are accepted. Thus, the total acceptance 
rate depends almost only on the Metropolis acceptance rate. We should note that a 
small core, such as the core of the A^ = 1442 surface under R/Ro = 10%, plays no 
role for the potential U{R). In fact, a new position X' on such small surfaces can 
move from X, which is outside of the core, to opposite side of the core. To check 
that a core is effective as U{R), it is helpful to see whether S i/N=\.5 holds or not 
at the transition, because the relation S i/N = 1.5 is violated only in the collapsed 
phase due to the potential U{R) as mentioned above. 

The total number of MC sweeps (MCS) after sufficiently large number of thermal- 
ization MCS is approximately 5 x 10^ ~ 7 x 10^ at the transition region of the 
surfaces of A^ = 8412, A^ = 10242, and N = 12252, and the total number of MCS is 
relatively small at the non-transition region and on the smaller sized surfaces. 



4 Results 

Snapshots of surfaces of N = 12252 are shown in Figs, ^a) and[21^b), which are 
separated by a first-order transition. Both surfaces appear to be in the smooth phase. 
However, the surfaces are different from each other because of the transition, which 
is the main result of this paper. The surfaces in Figs. [2a) and [2b) are respectively 
obtained at a = 760 and a = 730 under the condition R/Rq = 50%. The mean square 
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Fig. 2. Snapshots of the surface of size A'^= 12252 at (a) a -160 in the smooth phase and at 
(b) a = 730 in the collapsed phase. The ratio of R to Rq is R/Rq = 50%. The mean square size 
is about (a) = 175 and (b) = 145, which correspond to the mean radius Vl75 ^ 13.2 
and Vl45 - 12. The snapshots are drawn in the same scale. 

size is = 175 and X^ = 145 on the surfaces in Figs.[2ta) and[2tb), respectively, 
where X^ is defined by 

We should note that the values = 175 and = 145 correspond to the mean radii 
VT75 ^13.2 and VT45 ^ 12, both of which are relatively close to the core radius 
R = 10.26 inside the surfaces. 

The mean square size X^ in Eq. dH) is plotted against a in Fig. I2a), where R/Rq = 
50%. Figure Ob) is the variance of X^ defined by 

Cx^ = ^((^' - {X')f). (5) 

The peak values C™" vs. N are shown in Fig. [Sjc) in a log-log scale. The straight 
line in Fig. |3tc) is drawn by fitting the largest five data to 

^max ^ ^2/?^ 2J3 = 2.23 + 0.15. (6) 



Thus, we have ^C^f ~ A^-^ with = 1.17 + 0.08, and therefore, we find that the 
transition is of first-order because of the finite-size scaling theory [I36I37II . 

The Hausdorff dimension H is defined by 

X^ ^Ni (N ^ oo). (7) 

To obtain H at the transition point at R/Rq = 50%, we show X^ ™" and X^ ™' in 
Fig. Hd), where and X^'^°^ denote X^ in the smooth phase and in the col- 

lapsed phase, respectively. Since and x'^'^°^ are hardly obtained at the transi- 
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c 



H=2.26(2) 
collapse 



X2 



0.01 



0.005 



(b) |R/R,p0.5 



0:N= 


=12252 


A:N= 


=10242 


□ :N= 


=8412 


v:N= 


=4842 




600 



700 a 



X^/N 



0.02 



0.01 




max 
Cx2 



10" 



10" 



(c) 



R/Rn=0.5 



2(3=2.23(15) 



1000 



40 




. (f) R/i 


^0=0.3 

0— O N=4842 " 
N=2562 




r^N=1442" 



0.0004 



1/N 



800 



a 



900 



Fig. 3. (a) The mean square size X vs. a at R/Rq = 50%, (b) the variance C^j vs. a, 
(c) a log-log plot of C™ vs. N, (d) vs. N in the smooth and the collapsed phases at 
the transition point, (e) X^/N vs. and (f) X^ vs. a at /?//?o = 30%. The solid hnes 
connecting the symbols in (a), (b), and (f) are drawn as a guide to the eyes. 

tion point, we use the largest (smallest) at every N in Fig.[3ta) as X^'^™° (X^''^°'). 
By fitting the data excluding X^(A'^ = 6252) to Eq. (|7]), we have 



H = 2.18 + 0.02 (smooth), H = 2.26 + 0.02 (collapsed). 



(8) 



The value of H in the collapsed phase is very close to H in the smooth phase as 
expected. 

Although Si/N 1.5 due to the potential UiR), Z^ ^'/A^ is expected to be a non- 
zero finite constant in the limit of ^ oo. In fact, X^'^°^ is larger than the core 
radius squares and smaller than Recalling that both R^/N and X^'^'^°/N 

are independent of at ^ oo and that R^ < we expect that X^'™' is 

independent of at oo. Figure He) shows X^-^'/A^ and X^-'^^/A^ against 1/A^, 
where the straight lines are drawn by the linear fitting of the data excluding the data 
of A^=1442 and A^ = 6252. 



We see from Fig. [S^f) that X^ at R/Rq = 30% discontinuously changes even on the 
small surfaces. The transition point a^N) is not always identical with a where some 
physical quantities discontinuously change. A discontinuity seen X^ in Fig.[3];f) at 
R/Rq = 30% only implies a possibility that the surface configuration is trapped in a 
potential minimum. However, in this case the discontinuity clearly reflects a first- 
order nature of the transition. To the contrary, the transition point adN) is almost 
clear in Fig. 13 a) at R/Rq = 50% because physical quantities vary continuously, such 
as X^ of the A^^ = 12225 surface in Fig. ^a) at R/Rq = 50%, although a discontinuous 
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nature of the transition is unclear. 




sweeps 



4 xlO°0 
h(x') 



N=12252 
a=750 

J I I L 



2 4 6 xlO 
sweeps 



8 



0.5 



N=8412 
ot=720 




0.5 



(d) N=12252 
a=750 



100 110 120 130 




170 



Fig. 4. Variations of vs. MCS of (a) the = 8412 surface at a - 720 and (b) the = 1 2252 
surface at or = 750, (c) a normalized distribution /i(X^) corresponding to the variation in (a), 
and (d) h{X^) corresponding to the variation in (b). 



To see the discontinuity in more clearly on the surfaces at RjRo = 50%, we 
show variations of against MCS in Figs. Hl^a) andHfb). The discontinuity of X^ 
is hardly seen in the variation of X^ on the = 8412 surface, while it is clearly seen 
on the A^= 12252 surface. The corresponding normalized distribution h(X^) of X^ is 
shown in Figs.H^c) andHfd). A double peak structure can be seen in both of h{X^), 
and the peaks are clearly separated as increases. This implies a discontinuous 
transition. 




Fig. 5. (a) The bending energy S2INB vs. a at/?//?o-50% and , and (b) the variance £$2 vs. 
a, and (c) C^^^ vs. A'^ in a log-log scale. The symbol S2 is used to distinguish the extrinsic 
curvature energy from the intrinsic curvature one 5 2°'. 
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The bending energy 5* 2 is defined by 



S2 = J](l-nrnj), 



(9) 



where n, is a unit normal vector of the triangle i, and is the sum over all pairs of 
nearest neighbor triangles i and j. The surface shape is reflected in this 5 2, which 
is not included in the Hamiltonian. S2/NB vs. b is shown in Fig. [2a), where A^e is 
the total number of bonds. The variance of 5 2 is defined by 



(10) 



and it is plotted against a in Fig.|2b). The peaks C™" vs. are plotted in Fig. [3c) 
in a log-log scale. The straight line is drawn by fitting the largest five data to the 
relation 



/-■max 



2(r = 2.10 + 0.15. (11) 

Thus, we have cr= 1.05 + 0.08, where -^Jcf^ ~ N"'. This result implies a discontin- 
uous transition of surface fluctuations, and the result is consistent with that of Eq. 




2 4 xl0°0 

^^^^P^ h(S,) 



N=12252 
a=750 

J I I L 




2 4 
sweej3s 



6 xlO° 



1 



0.5 



N=12252 
a=750 




Fig. 6. The variation of 52/A^B vs. MCS of (a) the A^ = 8412 surface at a = 720 and (b) the 
N = 12252 surface at a = 750, (c) the normalized distribution h{S2) corresponding to (a), 
and (d) the normalized distribution h{S2) corresponding to (b). 

The variation of S2/NB against MCS is plotted in Figs.[6ta) and[6tb), where S2/NB 
reflects the shape fluctuations. A discontinuity of S2/NB is clearly seen in the vari- 
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ation of Fig. [6tb). Figures [6tc) and [6td) show the corresponding normalized his- 
tograms h{S2)- A double peak structure is clearly seen in h{S2) of Fig. Ilfd), while 
it is not in Fig. ^c). This also implies that the transition is of first order, because 
the double peak becomes apparent with increasing A^. 



S3/(3Nt) 



0.17 



0.165 



0.16 



N=8412 ^ 



h(S3) 




600 



700 



a 



0.16 0.17 
S4/(3Nt) 



Fig. 7. (a) The in-plane shear energy S3/{3Nt) vs. a a.tR/Ro = 50%, and (b) the normalized 
histogram h{S3) of S3/{3Nt) on the N= 12252 surface at a = 750. 



The in-plane shear energy 5 3 is defined by 



(12) 



where 0, is an internal angle / of a triangle, and Nt{= 2N-4) is the total number of 
triangles. 5 3 reflects an in-plane deformation of triangulated surfaces [i33J . and is 
not included in the Hamiltonian. 



FigureUta) shows S3/(3Nt) vs. a at R/Rq = 50%. The variation of 53/(3A^r) vs. a 
is also expected to be discontinuous at sufficiently large Adjust like the quantities S 2 
and X^, which characterize the out of plane deformation. The normalized histogram 
hiS^) is shown in Fig. |7tb), and we see a double peak structure in hiS^). This 
implies that the shape transformation transition can also be reflected in internal 
quantities of the surface. We note that the discontinuity is clearly seen in St,/(3Nt) 
at least when R/Rq < 30%. 

The intrinsic curvature energy 5^2^^ shown in Figs. [Ma)-Sc). As mentioned 
above, S™^/N does not always have a clear jump at the transition under i?/i?o = 50%. 
The reason why Sj"'/^^ appears to vary continuously against a under R/Rq = 50% 
is because the transition occurs many times during the simulations as we see in the 
variations of in Fig. Sfb) and S2/N in Fig. ^b). 

We should remark that S^"7A^ (A^ = 4842) at R/Rq = 50% in Fig. [Sla) is almost 
comparable to S^"7A^ (A^ = 4842) at R/Rq = 30%, R/Rq = 20%, and R/Rq = 0% in 
Figs.[8]^b)-[8]^c) in the collapsed phase close to the transition point. These values are 
pointed out by the left arrows in the figures. This observation allows us to consider 
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N= 12252 



s'";n 

0.0012 
0.00 
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a 
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Fig. 8. The intrinsic curvature energy S'f/N vs. a at (a) R/Rq = 50%, (b) R/Rq = 30% (c) 
R/Ro ^ 20%, 0%. The left arrows drawn at the ^ 4842 data indicate that S'^^/N in the 
collapsed phase is almost independent of the ratio R/Rq at the transition region. 

that the collapsed phase is degenerate. The collapsed states belong to a single value 
of S2^/N or to a narrow and continuous range of S^^N. This degeneracy comes 
from the fact that S is insensitive to the surface shape as we have discussed in 
Section |2l 

int 

So/N 




X :R/Ro=0 
A :R/Ro=0.2 
O :R/Ro=0.3 




X :R/Ro=0 
A :R/Ro=0.2 
O :R/Ro=0.3 



0.0004 0.0008 
1/N 



0.0004 0.0008 
1/N 



Fig. 9. S^2^/N vs. \/N at (a) the collapsed phase and (b) the smooth phase just below and 
above the transition points. The data are obtained on the surfaces of N = 4842, N = 2562, 
and N = 1442. The straight lines are drawn by the linear fitting of the data. In the limit of 
N ^ oo. 

smooth phase at the transition point. 



S2^/N is independent of R/Ro{< 0.3) in both of the collapsed phase and in the 



To see this degeneracy more clearly, we plot S^^N vs. 1/A^ in Fig. Ufa), where 
= 4842, = 2562, = 1442. The data 5^"7A^, including those just indicated 
by the left arrows in Figs. [Mb)-[S]^c), are obtained at the transition point. We see 
no dependence of S^2^/N on the conditions R/Rq = 0% ~ 30% and no difference 
in 5^7 even in the limit of N ^ oo. In the smooth phase, on the other hand, 
such a degeneracy is exactly expected in S^2^/N, because the smooth surface is not 
influenced by the core inside it. Figure |9tb) shows S^2^/N in the smooth phase. We 
should note that the size-effect is not expected in S2^/N. In fact, the size-effect can 
be seen only in the variances shown in Figs.[3tc) and[3c). 



Finally, we show S i/N vs. a in Figs. [TOlfa) - [TOl^ c). As mentioned in Section [2l 
the relation 5i/A^ = 3/2 can be violated due to the potential U{R). We see that 
S i/N discontinuously changes against a if R/Rq > 20% at least, while the relation 
5i/A^ = 3/2 holds at R/Rq = 0% and in the smooth phase at R/Rq < 30%. The 
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600 700 a 800 a 900 900 a looc 



Fig. 10. The Gaussian bond potential Si/N vs. a at (a) R/Rq = 50%, (b) R/Rq = 30%, (c) 
R/Ro = 20%,0%. 

variations of S i/N in Fig. \TO[ a) look very similar to those of S2/NB in Fig. I3a). If 
the gap in S i/N shown in Figs. ITOl h) and ITOl c) is very large, it can influence the 
surface size such as X^. However, the gap in S^/N is very small compared to the 
value of S i/N itself, and consequently, the influence of the gap in S JN on X^, or 
on the surface size, is negligible. Nevertheless, the discontinuous change in S i/N 
indicates that the surface area discontinuously changes. We should note that the 
gap inS i/N is seen on such small surface of N= 1442 at R/Rq = 20%. This implies 
that the core is effective as the potential U(R) in Eq. Q at least under R/Rq > 20% 
on the surfaces of N> 1442. Note also that the core slightly influences S i/N even 
in the smooth phase; in fact, S i/N at R/Rq = 20% in the smooth phase is slightly 
different from the one at R/Rq = 0% on the A^= 1442 surface. 



5 Summary and Conclusion 

In this paper, we numerically study a surface model of sphere topology focusing on 
whether or not the phase structure is influenced by a spherical core inside the sur- 
face. The model is an intrinsic curvature surface model, which is known to undergo 
a first-order transition between the smooth spherical phase and the collapsed phase 
on fixed-connectivity surfaces if it were not for the spherical core [i27il . 

We find that the model undergoes a discontinuous transition even when the surface 
encloses the core. The strength of the transition depends on the core radius R; the 
transition strengthens (weakens) with decreasing (increasing) R/Rq, where Rq is the 
radius of the equilibrium smooth surface at infinite curvature coefficient a ^ 00. 
The transition is of first-order when R/Rq < 50% at least. The results in this paper 
allow us to consider that the transition occurs at sufficiently large R/Ro(< 1) in the 
limit ofN^ 00. Moreover, we find also that the transition is non-trivially reflected 
in the Gaussian bond potential S i/N such that S JN discontinuously changes at the 
transition. The collapsing transition accompanies a discontinuous in-plane defor- 
mation such as the change of surface area. 
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We should comment on the reason why the transition at one R/Rq is identical to the 
transition at another i?/i?o- This is because the collapsed states are degenerate, and 
the degeneracy can be understood as follows: the transition looks dependent on the 
ratio R/Rq in the sense that the surface shape observed at one R/Rq is different from 
the one at another R/Rq. In fact, the collapsed phase, separated from the smooth 
phase by the transition, looks strongly dependent on the ratio R/Rq', the collapsed 
surface is almost indistinguishable from the smooth surface at R/Rq = 50%, while 
the collapsed surface is completely different from the smooth surface in the limit 
of R/Rq 0. However, the transitions are considered to be a single one, because 
all of those collapsed surfaces have almost the same value of the intrinsic curvature 
energy S2^/N. Thus, we see that the multitude of the transitions comes from the 
fact that S is insensitive to the surface shape. 

The results obtained in this paper imply that a large fluctuation of the density of 
membrane lipids is expected at the collapsing transition if the membrane encloses 
some artificial materials. In other words, the in-plane phase transition implies that 
the fluctuation of the lipid density, or of the distance between lipids, is divergent. 
Large fluctuations of these quantities are considered to be a signal of the phe- 
nomenon that the surface is almost broken. In fact, the obtained results indicate 
that a discontinuous change is observed in the area of surface, because S i/N is pro- 
portional to the area. Therefore, a discontinuous change of the density of vertices 
is also expected, because the surface area changes discontinuously while the total 
number of vertices is fixed. 

Finally, we comment on the future study. It is very interesting to study a surface 
model with a variable metric ^ on a two-dimensional surface M, from which a 
mapping X is defined such that X : M In this paper, we show that the 

phase transition is seen in an "intrinsic curvature" model, in which the surface area 
is not allowed to change in the whole region of the curvature coefficient a whenever 
the surface encloses no core. To the contrary, in such an "extrinsic curvature" model 
with a variable g, the surface area of X{M) is allowed to change at the collapsing 
transition even when the surface encloses no core. Therefore, the discontinuous 
change of S JN in the extrinsic curvature model with variable g is expected if the 
surface encloses a spherical core, and as a consequence, the in-plane deformation 
such as the change of surface area will be clarified in more detail. 
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